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Abstract The motion of a satellite around a planet can be studied by the Hill model, 
which is a modification of the restricted three body problem pertaining to motion of 
a satellite around a planet. Although the dynamics of the circular Hill model have 
been extensively studied in the literature, only few results about the dynamics of the 
elliptic model were known up to now, namely the equations of motion and few unstable 
families of periodic orbits. In the present study we extend these results by computing 
a large set of families of periodic orbits and their linear stability and classify them 
according to their resonance condition. Although most of them are unstable, we were 
able to find a considerable number of stable ones. By computing appropriate maps of 
dynamical stability, we study the effect of the planetary eccentricity on the stability of 
satellite orbits. We see that, even for large values of the planetary eccentricity, regular 
orbits can be found in the vicinity of stable periodic orbits. The majority of irregular 
orbits are escape orbits. 



1 Introduction 

The motion of a satellite around a planet, which revolves on a circular orbit around 

the Sun, can be modelled by the Hill's model of the three-body problem. Starting 

from the circular restricted three body problem (CRTBP), the Hill's approximation is 

achieved by translating the origin of the rotating reference frame to the planet and the 

1 

unit of length is scaled by the factor ji ' , where is the mass parameter of CRTBP. 
Then we let — > 0. The present study is restricted to the planar motion. Under these 
assumptions we obtain the Hill's equations (Szebehely, 1967) 

'i-2i^ = ii^ip-\ ^ + 2i = -7^p-^ (1) 

where p = (^^ + 77^)^/^. These equations admit an integral of motion, the well known 
Jacobi integral 

Cff = 3?V2p-i-(c2 + 772), (2) 
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Fig. 1 Families of periodic orbits of the Circular Hill problem. Thick (blue) or thin (red) 
curves indicate stable or unstable families, respectively. Grey regions indicate forbidden areas. 



and are invariant under the symmetries 



Although the circular Hill problem (CH) has only two degrees of freedom and is, 
moreover, autonomous, conservative and parameter free, it is not integrable (Meletlidou 
et al., 2001; Morales- Ruiz et al., 2005) and shows rich dynamics, which are depicted by 
computing Poincare surfaces of section (Henon, 1970; Chauvineau and Mignard, 1991). 

A detailed study of the CH problem has been given by Henon (1969,1970). He 
found the main families of periodic orbits, explored the phase space by using Poincare 
sections and computed the width of the stability regions. As we know, the families 
of periodic orbits consist the backbone of the phase space. The most important of 
them are presented in Fig. [T] (see also Henon, 2003; Batkhin and Batkhina, 2009). 
The domain ^ > corresponds to prograde and ^ < to retrograde motion. The two 
equilibrium (Lagrange) points Li and L2 are found for Cjj = 3'^^^ at the locations 
^ = ±3~^/^, rj^O. The distance Rh = 3"^^^ defines the HiU radius. We can see from 
Fig. E] that families of periodic orbits of retrograde satellites extend to distances much 
larger than the Hill radius and, interestingly, some of them are stable. 

The Hill's approximation can be applied in the same way to the spatial problem 
(Henon, 1974). Such a model has been used in space mission orbit design (Villac, 2003, 
2008). A single averaged model for the spatial Hill's problem has been also studied, 
providing particular solutions (see Vashkov'yak and Teslenko (2008) and references 
therein) . Hill models have also been derived for the case of binary satellites or asteroids 
(Henon and Petit, 1986; Chauvineau and Mignard, 1990). In these models the center 
of mass of the binary moves along a circular orbit around the Sun. A model where the 
center of the binary moves in an elliptic orbit has been given by Moons et al. (1988). 
This model uses a different formalism and aims to the study of the changes of orbital 
elements due to close encounters. An extensive list of references on the dynamics of 
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satellite motion is given in Waldvogel (1999). The general three body problem may be 
also used for studying satellite or binary motion (Hadjidemetriou and Voyatzis, 2011). 

The elliptic Hill (EH) model can be derived by considering the same Hill's approx- 
imation assumptions used in the CH model, but letting the planet move on an elliptic 
orbit around the Sun. Such a model has been introduced by Ichtiaroglou (1980, 1981) 
who also computed a few families of periodic orbits. A further study (Ichtiaroglou and 
Voyatzis, 1990) showed that all periodic orbits found were strongly unstable. It is ob- 
vious that the EH is a more appropriate model than the CH for studying the dynamics 
of a small satellite around a planet (or an asteroid) with eccentric motion. The aim of 
the present work is to study the dynamics of the EH problem in more detail, focusing 
on the main qualitative features of the phase space and orbit evolution. 

In the next section we describe briefly the EH model and discuss the existence and 
continuation of periodic orbits. In sectionOwe present the results of the computations 
of families of periodic orbits, giving particular attention to the stable ones. In section 
4 we explore the phase space of the model by computing various maps of stability 
and study the effect of the planetary eccentricity to the stability of the orbits. Our 
conclusions are given in section 5. 



2 The Elliptic Hill model 

Considering two primaries of masses mg (the Sun) and mi (the planet) that revolve 
on the plane Oxy in a Keplerian ellipse around their center of mass O, the motion of 
the planet along the rotating axis Ox with angular momentum Pg (see Hadjidemetriou 
1975) is described by the equation 

ii-Q^x:^'^ + {1- lifxr^"^ ^Q, (3) 

where ^ = mi/(mo 4- mi) and Q = (1 — n)Pg/jj,. By assuming the initial conditions 
a;i(0) = xio and ^(0) = 1, where 6 is the angular velocity, we get Q = Xiq. 

The motion of a massless body, which moves on the same plane Oxy, is described 
by the equations 

X — 2Qx^^y + Q'^x^'^x -— 2Qx^^ xiy -f ii{xx — x)r^^ — {jJ-xx + (1 — ^j)x) , , 

y = -2Qa;^^i -f Q'^x^'^y + 2Qx^^xix - ^lyr^^ - (1 - tJ-)yr2^, 

where r\ — {x — xi)^ + y^ and = [x + j^xi)^ + y"^ . Equations Q and (HJ are the 
equations of the restricted elliptic three body problem in the rotating frame. 
We apply the Hill's transformation 

1 i/3t 1/3 lr\ 

x = xi + ii'^, y = li'ri (5) 

to equations ^ and Q and, by letting ^ 0, we get (Ichtiaroglou, 1980) 

xi — xfgx^^ + x^'^ = (6) 



and 

^ = -2x{qX^ ^5 {-x^ " + x\qX^ * - p-") rj + 2xiQX^ -"xi^ 



i = 2xloXx 2»7 + (22:1 ^ + xi„x^ ^-p 3) 5 _ 2xlox^ 3^17? 
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with = + 77^. The equation (|6]) describes Keplerian motion with eccentricity, 
semi-major axis and period given by 



(1 -ep) y (1 -ep)3 

respectively. For bounded motion, xiq should be restricted in such a way that — 1 < 
Cp < 1. For ep > (ep < 0) the planet is at periapsis (apoapsis) at t = and it has 
ii(0) — 0. For Ep — {xi = xiq = 1) we obtain the equations HI of the CH model. 

Fixing the value of the planetary eccentricity Cp 7^ 0, the equations ([7}, which 
describe the motion of the massless satellite, constitute a periodic non-autonomous 
system of period T. The system ((Tjl possesses only the symmetry S and an orbit of 
initial conditions 

^(0) = v{o) = m = 0, m = m 

is a symmetric periodic orbit of period T' if it satisfies the periodicity conditions 

i{T'/2- Co, %) = 0, rj{T'/2; ^o, Vo) = 0. 

The period T' has to be an integer multiple of the planetary period T, i.e. T' — kT, 
1,2,... 

By varying the planetary eccentricity ep, a symmetric periodic orbit continues to 
exist, so that we get a monoparametric family (Ichtiaroglou, 1981). All periodic orbits 
that belong to such a family are isolated in phase space and have a period T' which 
depends on ep. For ep = 0, a family of the EH problem crosses a family of the CH 
problem. The crossing point must be a periodic orbit of the CH problem with period 
Tc = T' /X, where A = 1, 2, .. is the multiplicity of the periodic orbit. Since for ep = 
it is T = 27r, we conclude that the periodic orbits of the circular problem with period 
Tc — 2K7r/A can be considered as "bifurcation" points for the EH problem. Therefore 
these orbits are continued for Cp 7^ with multiplicity A and period T'=\Tc=2ktt. The 
ratio k/A is the resonance of the periodic orbit. 

Obviously, the number of bifurcation points on each family of the CH problem is 
infinite. In this study we restrict our attention to families that are of relatively small 
multiplicity and bifurcate mainly from stable periodic orbits. The linear stability of 
the orbits is determined by the two conjugate pairs of eigenvalues of the monodromy 
matrix of the system An orbit is stable (s) only when all eigenvalues lie on the 
unit circle. If one or two pairs of eigenvalues lie on the real axis we have single (u) or 
double (uu) instability, respectively. Finally when all eigenvalues are outside the unit 
circle (but they are not real) we have complex instability (cu). The different kinds of 
stability can be determined by computing stability indices as in Broucke (1969) (see, 
also, Ichtiaroglou and Voyatzis, 1990). A necessary (but not sufficient) condition for 
a family of the EH problem to emanate (from ep = 0) with a stable branch is the 
bifurcation point to be a stable periodic orbit. Certainly, the stability type can change 
along the family. 

In Fig. [2] we present the bifurcation points that are studied in the present work. 
The resonance k/A is indicated for each point. Most bifurcation points of Fig. [2] belong 
to the family / of retrograde orbits (left panel) and to the stable parts of families 
g, g of prograde orbits (right panel). The bifurcation points of families (;3, Kg and 
Hm play a special role, which will be discussed in the following section. We note that, 
as the Jacobi constant increases, the main families / and g tend to = and the 
Hill approximation degenerates into the two-body problem. In this case the bifurcation 
points to the EH model are of larger and larger multiplicity as ^ 0. 
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(a) (b) 

Fig. 2 Bifurcation points for continuation from the CH problem to the EH problem a from 
retrograde families b from prograde families of periodic orbits. The corresponding resonances 
are indicated. 



3 Families of periodic orbits 

A monoparametric family of the EH problem, which emanates from a k/X resonant 

At/A 

bifurcation point, will be denoted as f e , where F indicates the name of the family 
of the CH model where the bifurcation point belongs to. For the circular problem we 
adopt the naming introduced by Henon. The families a],^^, g'\!^ and /l^^ have been 
computed by Ichtiaroglou (1981) and are unstable (Ichtiaroglou and Voyatzis, 1990). 

Each orbit of an EH family is represented by a point in the three dimensional space 
^3 = {fip, ^Oi %} a-nd a family forms a characteristic curve in this space. In the following 
we will present, for convenience, the families in the projection plane (CoiEp)- Thick 
(blue) or thin (red) characteristic curves indicate stable or unstable orbits, respectively. 
The periodic orbits have been computed by the method of differential corrections and 
satisfy the periodicity condition with an accuracy of 10^^^. We always start from 
e-p — 0, where the period and the initial conditions are known from the CH model. 



3.1 Families /e - retrograte orbits 
3.1.1 The case k = 1 

We have computed the families up to multiplicity A = 9 and present them in Fig. [S] 
1/2 1/3 

The families fe and fe are unstable and extend outside of the Hill radius. The 
families with 6 < A < 9 start (at Cp = 0) as unstable but they posses a stable segment 
for moderate eccentricity values in both cases of the initial planetary position, namely 
for periapsis (ep > 0) and apoapsis {cp < 0). 

Families /^^^ and fe^^ in Fig. [3]seem to terminate abruptly at some critical points 
at Bp — ±0.215 and ep — ±0.185, respectively. However, as it is shown in Fig. [4l 
at these critical points, denoted as SN, saddle-node bifurcations take place. In these 
bifurcations the 1/1 resonant families Hgl^^^ and Hm^^^ , which bifurcate from the 
CH families Hg of multiplicity A = 4 and Hm of multiplicity A = 5, respectively, are 
involved. Particularly, in Fig. 2^ we obtain that the upper branch (cp > 0) of the family 
fe^'^ and the Hg\^^ emanate from the point SNi, while the lower branch (ep < 0) of 
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Fig. 3 Families /e'^^ of periodic orbits of the EH problem. Thick (blue) or thin (red) curves 
indicate stable or unstable orbits, respectively. 
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Fig. 4 a The continuation of family /e''* and its connection with families Hg\^^. b The 
continuation of family /i'^^ and its connection with family Hm\^^ . A closed path is formed. 
The symbol (uu) indicates parts of the families with double unstable orbits. 



the family /e and the HgJ emanate from SN2- In Fig. we see that the families 
/i^^ and Hml.^^ emanate from both critical points and SN/^, appearing for the 
eccentricity values ep — ±0.185. The two families connect smoothly and form a closed 
characterisric curve. 



3.1.2 The case k. = 2 

The computed families are presented in Fig. [5^. The families with A > 5 are unstable 
for the periapsis case but they start as stable when the planet is initially at apoapsis. 
The stable segment extends up to high eccentricity values in all cases. In contrast, the 



(a) 



(b) 



Fig. 5 a Families /e of periodic orbits of the EH problem, b The continuation of family 
fe^^ and its connection with family gSe^^ . Stability is indicated as (s) for stable (u) for unstable 
(uu) for double unstable and (cu) for complex unstable orbits. 




2 /3 

family /e is unstable for ep < and starts as stable for ep > 0. Then it seems to 
terminate but, as in the case of the family fe^^ , it connects with the family gSe^^ ■ The 
overall characteristic curve is given in Fig. [Sja. It forms a closed path, along which we 
find all possible types of linear stability. We note that the stable part is located outside 
the Hill sphere. 

3.1.3 The case k = 3 

We computed the families /e /e^^ and /e^^ presented in Fig.[6l All the corresponding 

bifurcation points are located outside the Hill sphere. The family /e^^ is unstable for 

small and moderate values of the planetary eccentricity, but it posses stable segments 

for \ep\ > 0.5. We note that for ep < the family becomes stable when it enters the 

3/5 

Hill sphere. On the other hand, the family /g is stable approximately in the interval 

—0.25 < ep < 0.25. Now, the family becomes unstable when it enters the Hill sphere 
3 /7 

for ep < 0. The family /g is unstable in all the explored domain. 
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3.2 Families ffe and g'g - prograte orbits 

In this case all bifurcation points and families are located inside the Hill sphere. 




5o 5o 

(a) (b) 

Fig. 7 a Families g'^^^ b Families g"^^^ ■ Thick (blue) or thin (red) solid curves indicate stable 
or simply unstable orbits, respectively, and dashed curves correspond to complex unstable 
orbits. 



3.2.1 Families ge' 

We computed the families with 5 < A < 9 and the families g^^^ with A — 11, 13 
and 15. They are presented in Fig. [7^. For ep > all families are unstable. For ep < 
the families g],^^ are also unstable but the families 3e^^ start us stable and become 
unstable when the planetary eccentricity exceeds a critical value, which is different for 
each case. 



3.2.2 Families g'^J 

The CH family g' has two branches (see Fig. [T} and each resonant bifurcation point 
appears twice, one in each branch; we use the notation "+" and " to distinguish 
between them. In this case the convergence of our numerical approach becomes delicate. 
The families we were able to compute are presented in Fig.fTJa. Also in this case we see 
that segments of stable orbits exist for low eccentricities. 



4 Phase space numerical exploration 

In this section we study the evolution of generally non-periodic orbits, associated with 
particular grids of initial conditions and various values of the planetary eccentricity, 
ep. The main orbit classification aims at distinguishing the orbits between bounded 
and escaping ones. We set the escape criterion to p > 15 Rh- Additionally, along 
the numerical integration of the orbit, we compute a Fast Lyapunov Indicator (FLI, 
Froeschle and Lega, 2000), in order to classify bounded orbits as regular or chaotic. 



10000 20000 30000 

time (t.u.) 

Fig. 8 The evolution of the FLI (without including the supremum) for some typical cases. 
Case (1) corresponds to a regular orbit, case (2) to a chaotic orbit and case (4) to a chaotic 
orbit that escapes for t < tmax- Case (3) refers to an unstable periodic orbit (see the text). 

The FLI, as defined by Voyatzis (2008), is given by tlie relation 

FLI= sup i|5(t)|, (9) 

wiiere S{t) is the deviation vector, whose evolution is given by the solution of the 
variational equations. In Fig. [S] we present the behaviour of the FLI for some typical 
trajectories. Case (1) corresponds to a regular (quasi-periodic) orbit, while case (4) 
corresponds to an irregular orbit that escapes at some t < tmax- Cases (2) and (3) 
are rather rare. Case (2) corresponds to a chaotic orbit which is sticky and does not 
escape for t < tmax- Case (3) corresponds to a trajectory that lies very close to an 
unstable periodic orbit. It is a typical case for all unstable periodic orbits presented 
in the previous section with relatively small planetary eccentricity (ep ^ 0.2). When 
the orbit follows the stable (unstable) manifold of the periodic orbit, the FLI decreases 
(increases) (see also Skokos et al., 2007). The overall (average) FLI evolution shows a 
linear increase in logarithmic scale, which is an indication of the existing instability. 

We have computed stability maps corresponding to a A'' x M grid of initial con- 
ditions. For each orbit numerical integration is performed up to tmax = 30000 time 
units or FLI < FLImax{=W^'^) or until the escape condition is fulfilled. Strongly 
chaotic orbits escape in relatively short time intervals. We observed that even if an 
orbit reaches the maximum FLI value, it shows a weakly chaotic evolution, at least for 

t <^ tmax - 



4.1 Henon stability maps 

In the circular Hill model, which has two degrees of freedom, the method of Poincare 
section can depict clearly the topology of phase space and the qualitative characteristics 
of the trajectories (Henon, 1970; Chauvineau and Mignard, 1991). For Ch > 3^^^ 
the zero velocity curves define closed regions on the 0^77 plane and all orbits, either 
regular or chaotic, are bounded. For C/f < 

34/3 

the majority of chaotic orbits are 
escape orbits. Henon (1970) determined the borders of bounded motion by considering 
the plane Ch — ^0 of initial conditions, with r;o = 0, and rja > 0, which are defined 
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C>1 



Fig. 9 Henon stability map for the CH model (ep = 0). The orbits of the grid are classified 
according to the color bar on the right of the map. The main families of periodic orbits are 
also presented with solid (stable) and dashed (unstable) curves. The two horizontal lines Rh 
determine the boundaries of the Hill's sphere. For $o = we get collision orbits (c.o.) while 
the grey regions (f.a.) indicate forbitten regions. 

by the Jacobi integral This diagram of bounded motion, which is called Henon 
diagram, has been used also by Shen and Tremaine (2008) in order to obtain bounded 
satellite motion by using numerical integrations of a 6-body model, which includes the 
Sun, the four giant planets and the satellite as a massless body. The Henon diagrams 
showed that regions of bounded motion are located near the main families (/, g and 
g') of stable periodic orbits. 

Following our methodology, the Henon map (or diagram) for the CH model (ep = 0) 
is given in Fig.[9l Empty (white) regions correspond to escape orbits. The regular orbits 
(dark regions) are distributed around the stable families of periodic orbits and this is 
in agreement with the results of Henon (1970). The reader can observe the tangle of 
stability around the stable family Hn, and this fact indicates the efficiency of the FLI 
method to provide detailed stability maps. Chaotic orbits (say orbits with FLI ^ 10^) 
are almost absent in the map, since the majority of chaotic orbits escape during the 
numerical integration. 

Using the same initial conditions as those of Fig. we set ep 7^ and repeat the 
computations. Note that positive or negative values of ep indicate that the planet is 
initially at periapsis (p-case) or apoapsis (a-case), respectively. Also, we remark that 
the points of the characteristic curves /, g and g' do not correspond to periodic orbits 
in this case. In Figs. [TUk and llOb we present the stability maps for the p and a case, 
respectively, and for Bp — 0.048 (Jupiter's eccentricity). We can observe that the width 
of the stable region of retrograde motion (i.e. around the family /) has shrunk but 
not significantly. The tail of this region {Ch < 0) seems now to be located below the 
characteristic curve / for the p-case and above it for the a-case. A rather significant 
effect of the non-zero planetary eccentricity value is observed in case of prograde orbits 
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Fig. 10 Henon stability maps of grid size 150 X 100 for the indicated planetary eccentricity 
values. Panels on the left correspond to the periapsis case and panels on the right to the 
apoapsis case. The color map is as in Fig. |9] 
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and mainly those located close to the g' family, where chaotic orbits appear. This is 
more clear in the a-case, where now the ordered orbits are obtained for higher values 
of Ch, compared to those of the circular case. 

By increasing the planetary eccentricity to Cp = 0.2 (periapsis case), we see from the 
corresponding map of Fig. llOfc that the area of the stability regions remains practically 
unchanged for retrograde as well as prograde orbits. However, for the apoapsis case 
(ep = —0.2) the stability region of retrograde orbits has shrunk significantly, while the 
stability region of prograde orbits has disappeared completely from the map (actually 
stability can be found for higher values of C^j). Finally in Figs. llOfe and II Of the case 
of a relatively high eccentricity value (|ep| = 0.5) is shown. It is remarkable that in 
the p-case we can still find regions of stable orbits for both retrograde and prograde 
cases. But retrograde stable orbits are found almost only inside the Hill's sphere. In 
contrast, the a-case shows an extensive reduction of the area of stable orbits. The only 
significant region of stable (retrograde) orbits left is located outside the Hill sphere. 

4.2 Stability maps along CH families 

The Henon maps reveal that, in the case of the EH model, regions of stable satellite or- 
bits are located around the stable families of the CH model. The effect of the planetary 
eccentricity to the periodic orbits of a given family, say A, of the circular model can 
be studied by constructing stability maps on the plane ~ Cpi where for each value 
we complete the set of initial conditions with % = Co = and 770 has the value that 
corresponds to the particular periodic orbit of the given family A. Thus, the initial 
conditions of the line Cp = correspond exactly to the family A of periodic orbits. 

The stability map along the family / is shown in Fig. 1111 Inside the Hill's sphere 
(— 3~^/^ < Co < 0) the stable orbits dominate and exist up to relatively high eccentric- 
ity values. Outside the Hill's sphere, we obtain a region of stability around the family / 
(ep = 0). Two of the gaps in this region are associated with the crossing of the unstable 
family g?> by the family / and can be also seen in the Henon maps of Figs. l^and llOb .b. 
This region is restricted to low (absolute) eccentricity values, because we have seen 
that, as the eccentricity increases, the stable regions located outside the Hill's sphere 
move away from the characteristic curve /. 

The stability map along the family g is shown in Fig. 112b ,. The point B indicates 
the pitchfork bifurcation of the family g' from the family g (Henon, 1969). No stable 
orbits are found around the unstable part of the family (on the right of the bifurcation 
point B). Around the stable part of the family, regions of ordered motion dominate, 
but the asymmetry between the periapsis and apoapsis cases is clearly seen. A similar 
asymmetry is also seen in the map along the family g', shown in Fig. 112b . So we conclude 
that the periapsis case gives significantly more stable prograde orbits, compared to the 
apoapsis case, in agreement with the results obtained from the Henon maps. 

4.3 Stability maps along EH families 

In the previous cases we associated regions of stable motion of the EH model with the 
families of periodic orbits of the CH model. Additionally, new regions of stability are 
expected to exist around stable periodic orbits of the EH model. We construct stability 
maps centered along a given family F of periodic orbits and defined by grids of initial 
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conditions on the plane ep-A(^Q. Each point of such a grid corresponds to the initial 
conditions + '^Coi % = Co = and t^q, where {^q, 7)q) are the initial conditions 
of the periodic orbit of family F at the planetary eccentricity ep, which has a period 
T' . Thus, the family F is located on the axis Z\^o ~ 0. For these maps we have set 
tmax = 20000T' and FLImax = 10^°. 

In Figs. [T5k .b we show the stability maps for the families /l^^ and /e^ '', respec- 
tively. Family /l^^ is stable in the segments —0.35 < ep < —0.23 and 0.23 < ep < 0.35. 
The stability map shows the existence of two islands of stability (II and 12) around 
these segments. For |ep| > 0.35, where the family is unstable, strongly irregular orbits 
exist, which escape after a few periods. In the interval —0.23 < ep < 0.23, where the 
family is also unstable, the existing instability gives a very thin chaotic zone in phase 
space, which is not clearly seen due to the resolution of the map. The wide stability 
regions (Rl and R2) that exist close to the unstable zone are associated with the stable 
family / of the circular problem, as we have mentioned above. 

2/5 

A similar dynamical situation is seen in the map of the family /e . This family 
is stable for —0.56 < ep < and close to this segment we can observe the existence 
of regimes (11-13) with stable orbits. For ep > the family is unstable, but for small 
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Fig. 13 Stability maps of grid size 150 X 100 along a. family b. family fl^^ . The color 
map is indicated in panel c . 

eccentricities the instability zone is very thin. In this instability zone the orbits are 
weakly chaotic and their FLI evolution is like that of Fig. |8] (case 3). The neighbouring 
stability regions Rl and R2 extend up to Cp ~ 0.25; for larger planetary eccentricities 
the orbits escape. 



5 Conclusions 

The elliptic Hill's model is an extension of the famous circular Hill's model. It is 
obtained by applying the classical Hill's transformation to the equations of the elliptic 
restricted three body problem in a rotating frame. In this paper we studied the main 
qualitative features of its dynamics by computing (i) families of periodic orbits and 
their stability and (ii) stability maps for the determination of regions with regular 
orbits. 

The families of periodic orbits in the EH model bifurcate from the families of the 
CH model and continue with parameter the planetary eccentricity, ep. We can define 
an infinite set of bifurcation points on the CH families (ep = 0) by increasing the 
multiplicity of the periodic orbits. We have determined a large set of bifurcation points 
that correspond mainly to stable periodic orbits, since we are interested in locating 
stable families. For each bifurcation point we computed the family of periodic orbits 
by varying the planetary eccentricity, ep. We classified the families according to their 
resonance and discussed their structure. We found that most of the families extend up to 
high absolute ep values for both the periapsis and the apoapsis case. However, in some 
cases the EH families reach a maximum value of |ep|, where a saddle- node bifurcation 
takes place. Then a new EH family is generated, which also meet a bifurcation point 
on a CH family. Nevertheless, there may exist EH families that do not originate from 
critical points of the circular problem, but such families cannot be computed by the 
approximation followed in this paper. 
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Most of the families are found to start as unstable, but in many cases family seg- 
ments exist having stable periodic orbits even for relatively high absolute eccentricity 
values (cp > 0.5). Generally, continuing the families for Cp > and Cp < we find that 

3/5 

at least one case gives simply unstable orbits. An exception is the family , which 
in both cases starts as stable. Family segments of double and complex instability also 
exist. 

In order to determine the phase space regions of regular orbits, we computed maps 
of dynamical stability defined on various planes. The majority of initial conditions cor- 
responds to orbits that are irregular-chaotic and finally escape. However, the Henon 
stability maps showed that around the stable families / and g of the CH model there 
exist wide regions of regular motion. As the planetary eccentricity increases, the sta- 
bility regions are affected significantly only when the planet is at apoapsis for t = 0. 
Prograde orbits are affected more than retrograde ones. 

The stable periodic orbits of the EH model are located at the "centers" of phase 
space regions with regular orbits. Instead, unstable periodic orbits are associated with 
chaotic orbits in their vicinity. For large absolute eccentricity values such chaotic orbits 
escape after a few planetary revolutions. For \ep\ ^ 0.2 the chaotic regions around the 
unstable periodic orbits are thin and are surrounded by regions of regular orbits. In 
these regions the orbits are trapped and do not escape, at least for relatively long time 
spans. 
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